WASP-19b: the shortest period transiting exoplanet yet discovered 

L. Hebb 1 , A. Collier-Cameron 1 , A.H.M.J. Triaud 2 , T.A. Lister 3 , B. Smalley 4 , P.F.L. Maxted 4 , 
C. Hellier 4 , D.R. Anderson 4 , D. Pollacco 5 , M.Gillon 2 - 6 , D. Queloz 2 , R.G. West 7 , S.Bentley 4 , 
B. Enoch 1 , C.A. Haswell 8 , K. Home 1 , M. Mayor 2 , F. Pepe 2 , D. Segransan 2 , I. Skillen 9 , S. Udry 2 , 

and 

P.J. Wheatley 10 
ABSTRACT 

We report on the discovery of a new extremely short period transiting extra-solar 
planet, WASP-19b. The planet has mass, M pl = 1.15 ± 0.08 Mj, radius, R pi = 1.31 ± 
0.06 Rj, and orbital period, P = 0.7888399 ± 0.0000008 days. Through spectroscopic 
analysis, we determine the host star to be a slightly super-solar metallicity ([M/H] = 
0.1 ± 0.1 dex) G-dwarf with T c g- = 5500 ± 100 K. In addition, we detect periodic, 
sinusoidal flux variations in the light curve which are used to derive a rotation period 
for the star of P ro t = 10.5 ±0.2 days. The relatively short stellar rotation period suggests 
that either WASP-19 is somewhat young (~ 600 Myr old) or tidal interactions between 
the two bodies have caused the planet to spiral inward over its lifetime resulting in the 
spin-up of the star. Due to the detection of the rotation period, this system has the 
potential to place strong constraints on the stellar tidal quality factor, Q' s , if a more 
precise age is determined. 

Subject headings: stars: planetary systems - techniques: radial velocities - techniques: 
photometric 
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Introduction 



Since the unexpected discovery of the first 'hot Jupiter', 51 Peg b (jMavor fc Quelozl Il995l ). 
exoplanets with an exceptionally wide variety of properties have been detected which have dra- 
matically changed our understanding of planetary physics. In particular, through the discovery of 
various tran siting planets, we have learned that extra-sol ar planets can ha ve radii much larger than 
Jupiter (e.g. iHebb et al.ll2009i ) or densities much higher (ISato et al.ll2005r). Many, b ut not all, 'hot 
Jupiters' have temperature inve rsions in their atm ospheres (e.g. lKnutson et al.ll2008l ). and they can 
have very low optical albedos (IRowe et al.ll2008h. Despite t heir short periods, not all transiting 
exoplanets have been tidally circularized (jGillon et al.ll2009al ). and both rocky (e.g. CoRoT-Exo-7, 
P ~ 0.85 days) and gas giant (e.g. WASP-12b, P ~ 1.09 days) planets can exist in extremely short 
period orbits. Here, we report on the discovery of a new extreme transiting extra-solar planet with 
the shortest orbital period yet detected which is on the verge of spiraling into its host star. This 
transiting planet can not only inform us about the properties and evolution of close-in planets, but 
it also has the potential to provide information about the characteristics of its host star. 

In this paper, we first describe all the observations that were obtained to detect and analyse 
the transiting star-planet system (§2)). We describe the data analysis in $3] where we present the 
planet and its host star. Finally in §U we discuss the implications of the planet's short period and 
its future evolution. 



2. Observations 



2MASS J09534008-4539330 (hereafter WASP-19) is an apparently unremarkable 12th mag- 
nitude (V = 12.59), G8V star in the southern hemisphere located at a =09:53:40.08, 5 =- 
45:39:33.0 (J2000). The target was obser ved with the WASP-South telescope and instrumentation 
(jPollacco et al.1 120061 : 1 Wilson et al.1 120081 ) in the winter and spring observing seasons from 2006 to 
2008. 1496 photometric data points were obtained between 4 May - 20 June 2006, 6695 measure- 
ments were made between 18 Dec 2006 - 18 May 2007, and 8968 observations were taken from 18 
Dec 2007 - 22 May 2008. All data sets were p rocessed independently with the standard WASP 
data reduction pipeline and photometry package (jCollier Cameron et al.ll2006l ) . The individual data 
points have typical uncertainties of ~ 0.02 mags including poisson noise and systematic noise. The 
resulting light curve s were then run through our implementation of the box least squares algorithm 
(jKovacs et al.ll2002l ) designed to detect periodic transit-shaped dips in brightness. 



The target was initially flagged as a transiting planet candidate because a strong periodic 
signal was detected in the the 2007 data. The phase-folded light curve showed a square-shaped dip 
in brightness with a depth, 5 ~ 25 mmag and duration, r ~ 1.2 hours, consistent with a planet 
sized object around a main sequence star. Further, a periodic transit was also apparent in the 
2006 data when phase-folded with the 2007 ephemeris, and a transit was subsequently detected in 
the 2008 season of data. Therefore, we classified the object as needing follow up photometry and 
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spectroscopy to assess the planetary nature of the system. The phase-folded light curve containing 
all WASP-South data is shown in Figure [TJ 

WASP-19 was observed photometrically with the 2m Faulkes Telescope South (FTS) on 17 
December 2008 during transit. 139 Pan-STARRS z-bancfl observations were made over 3.3 hours. 
The images were observed in 2 x 2 binning mode such that one binned pixel corresponds to 0.279". 
They were processed in the standard way with IRAF using a stacked bias image, dark frame, and 
sky flat. Minimal fringing was present in the z-band images due to the deep depleti on CCD in the 
camera, so no fringe correction was applied. The DAOPHOT photometry package (jStetsonlll987l ) 
was used to perform object detection and aperture photometry with an aperture size of 8 binned- 
pixels in radius. The 5' x 5' field-of-view of the instrument contained 53 comparison stars that 
were used in deriving t he differential m agnitudes with a photometric precision of 1.3 mmag. We 
measured the red noise (jPont et al.ll2006l ) in the light curve on a 30 minute timescale to be 449 ppm 
and added this value in quadrature to the formal uncertainties on each data point. The resulting 
light curve is shown in Figure [2j 

Thirty- four radial vel ocity measurements were obtained with the CORALIE spectrograph on 
the 1.2m Euler telescope ([Baranne et al.lll996l ; IWilson et al.ll2008l ). The stable, temperature con- 
trolled, high-resolution echelle spectrograph has a resolution of R ~ 55000 over the spectral region 
from 3800-6800 A. The WASP-19 spectra, obtained between 29 May 2008 and 23 April 2009, were 
processed through a slightly updated v ersion of the CORALI E data reduction pipeline. In addition 
to the standard pipeline described in iBaranne et al.1 (jl996l ). we corrected for the blaze function 
and scaled the fitted cross-correlation region to match the full-width half maximum of the object. 
The final radial velocity (RV) values were obtained by cross-correlating the spectra with a G2 tem- 
plate mask. Table [1] presents the radial velocity measurements of WASP-19 at each Bary c entric 
Julian date, the la Poiss on errors on the velocities, and the line bisector span measurements (jGrav 



1988 



Queloz et al.ll200ll ). Based on our experience, we adopt uncertainties on the line bisector 



measurements of twice the measured RV errors. 

The radial velocity of the star varies sinusoidally with the same period measured from the 
photometry (see Figure [3]). In addition, the line bisector spans, which are used to discriminate spot 
induced velocity variations and the effects of line-of-sight binarity, show no correlation with radial 
velocity within the uncertainties. The slope of the bisector versus RV (Figure 0]) is —0.006 ± 0.037, 
and the bisector measurements have an r.m.s. scatter of ~ 50 m s . 

Although all existing photometric and spectroscopic data suggests WASP-19 is orbited by a 
short period transiting extra-solar planet, we explore the possibility of a false positive detection. In 
general, it is difficult to mimic the photometric and spectroscopic observations of a transiting planet 
without showing a visible second star in the spectrum, a significant bisector trend, and/or incon- 
sistencies between the transit duration and host star spectral type. Starspots can cause periodic 



1 http:/ /pan-starrs. ifa.hawaii.edu/public/design-features/cameras. html 
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Table 1: Radial velocity measurements of WASP-19 obtained with the CORALIE spectrograph. 
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low-amplitude RV variations (e.g. iHuelamo et al.l 120081 ; iDesort et al.l 120071 ) but not photometric 



transits as well, thus a single star blended with a fainter stellar eclipsing binary (EB) is the pre- 
ferred false positive scenario for transiting planets. Unfortunately, no comprehensive simulations 
have been performed which model the expected RV variations, eclipse shapes, and bisector slopes 
for different blended EB scenarios, and performing such simulations is beyond the scope of this 
discovery paper. Instead, we explore the blended EB scenario through qualitative reasoning. 

In order to produce a flat-bottomed, 2% transit, as seen for WASP-19, the flux ratio of the 
EB compared to WASP-19 would have to be small enough that the EB was undetectable as a peak 
in the cross-correlation function. However, it could not be so small that the necessary unblended 
eclipse depth would require nearly equal sized EB components and therefore, a V-shaped eclipse (i.e. 
0.05 < Feb/ Fwi9 < 0.2). Although the eclipsing star would have to be small (and presumably 
less massive) compared to its primary to create the flat-bottomed eclipse, the RV amplitude of 
the visible EB primary would still be 40-90 km s _1 for all reasonable mass ratios due to the 
high inclincation angle needed to eclipse, the short orbital period of the observed transits, and 
the relatively deep transit (i.e. brown dwarf mass eclipsing objects would produce much shallower 
transit depths). Therefore, the RV variations of the EB could not be hidden within the WASP-19 
spectral features given the resolution of the CORALIE data. 

Furthermore, a short period stellar EB would almost ce rtainly be tidally sy nchronized with 



40 - 80 km s -1 . In the analysis of HD 41004 A (jSantos et all 120021 ). which exhibits 
planet-like RV variations due to a blended M dwarf + brown dwarf spectroscopic binary (SB), the 
bisector correlation is most dependent on the width of the visible SB component. This system, 
in which the SB is only 3% as bright as the single star and the width of the SB cross correlation 
function is ~ 8 km s -1 , shows a significant bisector correlation (slope of 0.67). According to their 
simulations, higher rotational broadening would produce an even greater bisector slope. Therefore, 
we would expect WASP-19 to show a significant bisector trend if the RV signal was caused by a 
rapidly rotating, blended EB, rather than a transiting planet. 

Finally, we see no difference in the eclipse depth of the z-band FTS transit compared the 
SuperWASP transit taken in a bluer, V+R filter. Although this is not a strong constraint given 
the scatter in the SuperWASP photometry, an eclipse by a stellar object could show eclipse depth 
variations in different filters which we do not see. In summary, we conclude that the existing 
photometric and radial velocity variations of WASP-19 are most likely due to the presence of a 
transiting extra-solar planet. 



3. Analysis 

A combined analysis of the radial velocity curve and light curve of a transiting planet host star 
will provide direct measurements of the mass and radius of its orbiting planet with one additional 
constraint, e.g. the stellar mass. Below, we describe the determination of the stellar mass and other 
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host star properties. 



3.1. Spectroscopic parameters 



The individual CORALIE spectra were co-added into a single higher signal-to- noise (S/N) 
spectrum which was then used to measure the stellar temperature, gravity, metallic ity, v sini and 



eleme ntal abundance information through comparisons with the synthetic spectra of lCastelli et al 



(jl997l ). The results of the analysis are listed in Tabled! The quoted uncertainties directly correlate 
with the modest signal-to-noise of the co-added spectrum (S/N ~ 70). 



The spectral synthesis technique is described in detail in iGillon et al.l ()2009al ). therefore we 
only give a short summary here. The Ha line was the primary temperature determinant, while 
the Na I D and Mg I b lines were used as surface gravity (log g) diagnostics. By measuring the 
equivalent width of several clean and unblended metal lines, we derived abundances for the elements 
listed in Table [2j Due to the spread in abundance measurements, we adopt an overall metallicity 
for the system of [M/H]= 0.1 ± 0.1. The microturbulence (1.1 ± 0.2 km s^ 1 ), which directly affects 
the abundance measurements, was derived from the Fe I lines using Magain's (1984) method. 

The projected rotational velocity (usini) was determined by fitting the profiles of several 
unblended Fe I lines. We accounted for line broadening due to the instrumental FWHM (O.llA) 
which was determined from tell uric lines ar ound 6300A, and the macroturbulence (2 km s _1 ) which 
was based on the tabulation bv lGravi (|2008l ). Finally, we do not detect Li I in the stacked spectrum, 
and can therefore only put an upper limit on the abundance of this element of log^4(Li) < 1. 



3.2. Host star mass 



By comparing the effective temperature, metallicity, and mean stellar density (p*) of WASP-19 
to theoretical stellar models, we determine its mass. The stellar density i s dependent on the shape 
of th e transit and largely independent of any assumptions or models (jSeager fe Mallen-Ornelas 
2003). Thus, after deriving the spectroscopic parameters, we model the WASP-South and FTS 
tr ansit light curves of WASP-19 using the Markov-chain Monte Carlo (MCMC) routined described 



m 



Collier Cameron et al.l (120071 ) to derive the stellar density and its uncertainty. The code uses 



the MCMC approach to simultaneously solve for the orbi tal and physical pr operties of the star- 
planet system. We apply the limb darkening coefficients of lClaretl (|2000i . 120041 ) for the appropriate 
temperature of the star and wavelength of the light curves. We note that the eccentricity value 
has an effect on the stellar density determination, thus we first allow the eccentricity to be a free 
parameter. Since the resulting eccentricity value gives a non-signficant 1.5a detection, we solve for 
the stellar density a second time while fixing the eccentricity to be zero. 



Interpolating among the 



Girardi et al 



(|2000l ) stellar evolution tracks as described in 



Hebb et al 
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Table 2: Stellar properties of WASP-19 obtained from the NOMAD catalogue and derived from our 
analysis of the spectra and light curves. 

Parameter WASP-19 



RA(J2000) 09:53:40.08 

Dec(J2000) -45:39:33.0 

J 10.911 ±0.026 

H 10.602 ± 0.022 

K 10.481 ± 0.023 

Ura —41.3 ± 2.5 mas yr~ 

/^DEC 16.5 ± 1.9 mas yr~ 

U -49t?5 km s_1 

1-2 
-2 
-7 
-10 



V -25l| km s" 1 

W -13+L km s- 1 



T cff 5500 ± 100 K 

log g 4.5 ± 0.2 

& 1.1 ± 0.1 km s" 1 

4 ± 2 km s- 1 



v smi 



[Fe/H] 0.02 ± 0.09 

[Si/H] 0.15 ± 0.07 

[Ca/H] 0.12 ± 0.15 

[Ti/H] 0.13 ± 0.12 

[Ni/H] 0.10 ± 0.08 

log A(Li) <1. 



1.13 ±0.12/90 
P rot 10.5 ± 0.2 days 
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( 2009 ). we derive a mass for the host star of M* = O.95±^M using the circular orbit solution. 
The non-circular orbit value for the stellar density results in a mass of 0.96 M@, well within the 
existing error bars (adopted from the non-circular solution) . Figure [5] shows a plot of the position 
of WASP- 19 in a modified Hertzprung- Russell diagram as compared to the theoretical tracks. The 
errors on the stellar mass are dominated by the uncertainty on the metallicity. Here, we take 
into account uncertainties on the stellar temperature, metallicity and density. We do not include 
systematic uncertainties on the chosen stellar model whi ch are difficu l t to d etermine, but we might 
expect this additional uncertainty to be at least ~ 4 %. ISouthworthl (120091 ) finds that with stellar 
parameters measured at the limit of our technological and theoretical ability for HD 209458b, the 
variation in the mass determination using four different stellar models is ~ 4%. 



3.3. Host star age 

The isochrone fitting also allows for deriving an age estimate for WASP-19. According to these 
stellar models, WASP-19 is a main-sequence star with an age of 5.5I45 Gyr when we adopt the zero 
eccentricity value for the stellar density. This essentially places a weak constraint on the stellar age 
to be > 1 Gyr. If we adopt the value for the stellar density from the solution when the eccentricity is 
a free parameter (non-significant 1.5a result), we are unable to place any constraints on the stellar 
age from the isochrones. The depletion of lithium in the atmosphere of a G8V star can also be 
used as an age indicator. However, the non-detection of Li 1 in WASP-19 agai n gives only a weak 



const raint on the age suggesting the star is older than the Hyades (0.6 Myr, ISestito Randich 



2005) 



We, therefore, examine the three dimensional velocity of WASP-19 as compared to theoretical 
Galactic model stars to estimate the probability that WASP-19 is part of a young disk population. 
Using the catalogue proper motions and measured systemic radial velocity, we calculate its U, V, 
and W space motion (given in Table [2]) compared to the Sun. We estimate the distance of WASP- 
19 to be 2 50+gjj pc using the measured spectral type, G8V, with absolute magnitude, M v = 5.6 



fromlGrayj (|1988l ) and V magnitude from the NOMAD catalogue (V = 12.59). We adopt a ±0.2 
uncertainty on the magnitude measurement as it is not given in the catalogue and a generous 
error of ±2 subclasses in spectral type to determine the uncertainty on the distance estimate. 
We select a set of model main sequence stars of spectral types F7-K7 in a small volume around 
WASP-19 (I = 27 3 ±3°, b = 7 ± 3°, distance= 150-400 pc) using the Besancon Galactic model 



([Robin et al.ll2003f ). We generate 100 resolutions of the simulation which provide the heliocentric 
velocities, metallicities, and population classes for over 500,000 model stars. Of the model stars 
with metallicities of 0.0-0.2 dex and with space motions within the errors of the calculated U,V, 
and W values of WASP-19, 35% of the model stars have ages of < 1 Gyr (population class 1 and 
2), 59% have ages of 1-5 Gyr (population class 3-5) and 6% have ages of 5-10 Gyr. Contributions 
from the thick disk, halo and bulge populations are negligible. This analysis suggests WASP-19 
has a 65% probability of being older than 1 Gyr. 
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In summary, we present three different age dating techniques which all suggest WASP- 19 is 
older than ~ 1 Gyr, however a precise age for the star cannot be determined from the existing 
data. 



3.4. Host star rotation period 

We search for variability in the WASP-South light curves of WASP-19 caused by asymmetric 
starspots on the photosphere which modulate the flux. Any starspots will rotate in and out of 
view with the stellar surface, such that the period of the variability gives the rotation period of 
the star. To exhibit rotational variability, the star must have a sufficient coverage of starspots to 
create a variable brightness signal which is detectable given the scatter in the photometric data. 
Furthermore, starspots evolve on timescales of weeks or months causing the amplitude and phase of 
the variability to change, therefore, each season of WASP-South data was examined independently. 

To detect the rotational variability, we determine the improvement in x 2 over a flat, non- 
variable model when a sine wave of the form y = ao + a\sin(ujt + 03) is fit to each season of the 
WASP-South data phase-folded at a set of trial periods, P ro t = 2tt/uj. We subtract all transits 
from the light curves using the model derived from the parameters in Table |3] before fitting the 
sine curve model. We test periods between 0.2 — 50 days and find a strong periodic signal in the 
2007 data with P rot = 10.5 days and amplitude, a\ = 7.6 mmag. The phase-folded light curve and 
periodogram of normalized Ax 2 values are shown in Figure |6] and Figure [7| (top) , respectively. In 
the 2008 data, we also detect a weaker signal with a similar period, P ro t = 10.6 days, and with an 
amplitude of 3.6 mmag. 

To assess the veracity of the sinusoidal signal detect ed in the 2007 data, we determ ine the sig- 



nificance and the false alarm probability (FAP) following I Zechmeister fc Kiirsterl |2009) (employing 
the residual variance normalization). The equations include, N, the number of independent data 
points, and M, the number of independe n t freq uencies, as well as the measured peak value in 



the periodogram. According to ICummingi (J2004J), the number of independent frequencies can be 



approximated by the duration of the time-series data times the difference between the highest and 
lowest frequencies tested (here M = 753). The number of independent data points depends on the 
level of red noise in the light curve. We empirically determine N by generating a new light curve 
in which the red noise is preserved, but the periodic signal is destroyed. We randomly re-order 
the individual nights of data, so that the integer part of the light curve time values are shuffled, 
but the fractional parts (containing the red noise) remain the same. We then run the sine fitting 
program on the shuffled light curve and find N = 664 by assuming the highest peak in the resulting 
periodogram (Figure [7| (bottom)) has a 95% probabi lity of being false (FAP = .95). Using our 



calculated N and M, we then apply the equations in IZechmeister fc Kiirsterl (|2009l ) to the results 
of the sine fitting on the original 2007 light curve and find a highly significant periodic signal with 
a Prob(p > Pbcst) = 1-4 x 10 -10 and FAP = 1 x 10 -7 . Thus, we adopt a rotation period for 
WASP-19 of P rot = 10.5 days. 
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We calculate the error on this measurement using the formula in iHorne Baliunasl (|1986l ) and 
find op = 0.08, but given the variation in period values we measure using two other techniques 
(Lomb-Scargle and auto-correlation), we find this to be underestimated and suggest ap = 0.2 days 
is a more realistic uncertainty. Finally, we note that the 10.5 day rotation period of the star 
measured via the photometry is consistent with the less precise fsini value of 4 ± 2 km s -1 which 
corresponds to a rotational period of between 8-24 days (assuming the spin axis of the star is 
perfectly aligned with the orbital axis). 

The ro t ation period of a main sequence star can also be used to estimate its age. According to 
the iBarnesI (|2007l ) gyrochronology relationship, a rotation period of P ro t = 10.5 days corresponds 
to an age of 500-600 Myr for a G8V star like WASP-19 with B-V=0.74. This value is not consistent 
with the older age inferred from the isochrones, lithium, and space motion. In § [H we discuss the 
implications of the possible discrepancy between the different age indicators. 



3.5. Planet properties 



We solved for the properties of the star-planet system by running an MCMC code (jCollier Cameron et al 



20071 ) using as inputs the WASP-South light curve, the FTS light curve, the CORALIE radial ve- 
locity curve, and the stellar mass. Initially, we allowed the ecosw and esinw to be free parameters, 
but the low, non-zero values that resulted were not significant. Therefore, we also solved for the 
parameters of the system while forcing the planet to be on a circular orbit. The difference in stellar 
and planet properties derived for the circular and non-circular orbit cases is negligible and within 
the la uncertainties on all parameters. However, for completeness, we provide the resulting solu- 
tions for both cases. Table O gives the final model parameters for the transit and radial velocity 
curves as well as the physical properties of the star-planet system. 

Note, we do not remove the rotational variability from the light curves before deriving the 
final system parameters. The model fit is dominated by the FTS data, and the low-amplitude 
10.5 day rotational variability is negligible on the 3.3 hour timescale of this light curve. We confirm 
this by fixing the orbital period to the value found when analysing all the light curve data and 
fitting just the FTS light curve and radial velocity data. All the resulting light curve parameters 
and physical system parameters are well within their reported la uncertainties. Therefore, we 
do not 'pre-whiten' the WASP-South data before performing the final transit model fits on all 
the existing photometric data. However, it is important to note that we cannot determine if any 
starspots were present on WASP-19 during the FTS observations which might affect the resulting 
parameter determinations because of the relatively small amount of out-of-transit data. Additional 
high quality transit data would allow for investigating any variations in the derived parameters due 
to starspots. 



F inally , to derive the final parameters, we use the 4-coefficient limb darkening model by [Claret 

(|2004l . l2000l ) with the ATLAS model atmospheres, an effective temperature of 5500 K, and a log go! 
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4.5. However, we tested a range of limb darkening coefficients with different temperatures, log g 
values, and model atmospheres and found the resulting parameters to be highly robust to changes in 
these coefficients. We suspect this is due to the fact the FTS z-band data dominates our parameter 
results, and this filter is not as susceptible to limb darkening uncertainties as bluer filter data. 



3.6. Transit Timing 



We measured the heliocentric Juili a n dat e of the mid-transit times for WASP- 19b using the 
technique described in I Anderson et al. (200 9J) to search for variations which could indicate the 
presence of a outer planet (e.g. lAgol et al.1 120051 ) . In general, the SuperWASP transits are not 
precise enough to provide a useful constraint on the existence of a third body in the system, but 
these data do rule out transit timing variations larger than ~ 15 minutes. 



Discussion 



The most striking aspect of WASP-19b is its extremely short orbital period. With a period, 
P = 0.7888399 ± 0.0000008 days, WASP19b is the shortest period planet yet discovered. What 
physical processes in the evolution of the planet have lead to such a close separation? Did the initial 
migration process leave the planet in its extremely short period orbit or has there been subsequent 
evolution of the orbital separation? The observations presented here suggest the possibility that 
WASP-19b has been spiralling into its host star throughout its lifetime and has spun up its host 
star in the process. 

Most transiting extrasolar planets will ultimately spiral into their host stars becau se there is 



insufficient total angular momentum i n the systems to reach a state of tidal equilibrium (jHut 



1980 



Rasio et al,lll996l : iJackson et al.ll2009l ). However, the timescale for this evolution is not well known 
because the stellar tidal quality factor, Q ' , is uncertain to at least three orders of magnitude. Values 
between 10 6 < Q' s < 10 9 are reasonable ( Jackson et al.ll2008l . and references therein). Furthermore, 
the planetary tidal quality factor, Q' is equally uncertain, but it affects the evolution of the stellar 
spin and planetary orbital separation to a much lesser degree. 

For WASP- 19b, the ratio of total to critical angular momentum is well below the limiting 
value of one, thus ensuring the planet will ultimately collide with its host star, but the lifetime of 
this evolution ranges from 4 Myr to 4 Gyr depending on the value used for Q' s . However, unlike 
most other transiting planets, we have measured the rotation period of WASP-19 which contributes 
additional information that can be used to place tighter constraints on the overall lifetime of the 
planet and on the tidal quality factor of the star. 



By integrating the coupled equations of orbital eccentricity and separation (jDobbs-Dixon et al. 

2004) while conserving the total angular momentum of the system, we can model the future evo- 
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Table 3: WASP-19 system parameters and la error limits derived from the MCMC analysis. 



Parameter 



Symbol Value (e fixed) 



Value (e free) Units 



Transit epoch (BJD) 
Orbital period 
Planet/star area ratio 
Transit duration 
Impact parameter 



T 
P 

(Rp/R s ) 2 

Pr 
b 



0.0001 
-0.0002 
-0.0000008 
-0.0000008 
0.0004 
0.0004 
0.0006 
0.0006 

0.62+r 3 



2454775.3372 
0.7888399 

0.0203; 

0.0642" 



+0.0002 
-0.0002 
-0.0000008 
-0.0000008 
0.0004 
0.0004 

0.0643+ - 0006 



2454775.3372 
0.7888399 
0.0203! 



-0.03 



-0.0007 

62+ 03 

u - DZ -0.03 



days 
days 

days 



Stellar reflex velocity 
Centre-of-mass velocity 
Orbital semimajor axis 
Orbital inclination 
Orbital eccentricity 
Longitude of periastron 
eccentricity x cos (a;) 
eccentricity x sin (a;) 



7 
a 
I 

e 

ecosco 
esinw 



20.785341858" 



0.0165 



-0.0005 
-0.0006 



8o.5i8:? 

(fixed) 
(fixed) 
(fixed) 
(fixed) 



0.256 



+0.005 
-0.005 



20.78535 
0.0164 



+0.0003 
-0.0003 
-0.0005 



0.0006 



80.8181 

0.0218:8? 

+112 



-76 
0.004 



-23 
+0.009 
-0.009 



-0 02 +a02 
u - uz -0.02 



km s 1 
km s -1 
AU 

degrees 
deg 



Stellar mass 
Stellar radius 
Stellar surface gravity 
Stellar density 



i?* 

log g* 
p* 



0.96 
0.94 
4.47 
1 



0.09 
0.10 
0.04 
0.04 
0.03 
-0.03 

1 O+0.09 

lo -0.09 



95 +0 - 10 
u - yo -o.io 

n OQ+0.05 
u - yo -0.04 
A zio+0.03 
4 - 48 -0.03 

1 1O+0.12 
i - ly -0.11 



M 

Rq 

[cgs] 

p® 



Planet radius 
Planet mass 

Planetary surface gravity 
Planet density 

Planet temperature (A = 0,F= 



R p 
M p 
logfi-p 
P P 
Teq 



1.31 
1 



0.06 
0.06 
i cr+0.08 
1 J_n i 



-0.08 
o iq+0.03 
°- ly -0.03 
H0.06 
-0.05 



0.511 



20091^ 



1 no+0.07 

1 14+0.07 
i - i4 -0.07 
o 9n +0.03 
°- zu -0.03 

u - O4 -0.06 

1993±|| 



Rj 

Mj 
[cgs] 

PJ 
K 
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lution of the planet as it spirals inward and further spins up the host star. In the analysis, we 
include magnetic braki ng which takes the form C l = —ntl 3 with k = 3.88 x 10~ 7 s as evaluated 



using equation (12) of ICollier Cameron &: Jiankd (|1994h . We use the current orbital separation 
derived from the MCMC analysis (Table [3|) and the measured rotation period {P TO t = 10-5 days) 
as initial conditions. Figure [8] shows the future evolution of the planet in orbital separation, and 
in Figure we plot the corresponding rotational evolution of the star. 

The results show that in the Q' s = 10 9 case, the tidal interaction is so weak that the star's 
future evolution is dominated by magnetic braking. In this situation, the relatively short rotation 
period could not have been caused by prior tidal interactions with the planet and would instead 
suggest a young stellar age similar to that of the Hyades (~ 600 Myr). Although it is possible 
WASP-19 is as young as the Hyades, the isochrones analysis, the non-detection of lithium, and the 
stellar velocity all favor an older age for the star. Therefore, a lower value for Q' s is preferred. 

For all models with Q' s < 10 8 , the star's spin period has already passed through a maximum 
and is decreasing as the planet spirals in. However, in the Q' s = 10 6 case, the remaining lifetime 
of the planet is extremely short (4 Myr). Therefore, the existing data appears to favor values of 
Q' s = 10 — 10 8 which are high enough to give significant life expectancy, but low enough to have 
reversed magnetic braking and initiated spiral-in. This scenario would allow for the extremely 
short period of the planet to be a consequence of further evolution in orbital separation after the 
initial migration process early in its formation and evolution. However, the Q' s would have to 
be 1-2 orders of magnitude greater than the nominal value of 10 6 that is typically adopted (e.g. 



Bodenheimer et al. 



2003 



Jackson et al.ll2008l ). 



It is important to note that Q' s is a simple parameterization of the complex physics involving 
the interaction between the tides raised by the planet on its star and the turbule nt viscosity in the 



stellar convection zone and inertial wave modes excited in the stellar interior (Rasio et al. 1996 



Sasselovil2003l ; lOgilvie fc Linll2007l ). A more detailed analysis of this system and others like it (e.g. 
WASP-18, OGLE-TR-56b, WASP-12) will hopefully lead to a better understanding of the physics 
modulating tidal dissipation in stars, since currently, there is no comprehensive theory which is 
able to explain observations of bot h main sequence binary stars and close-in extra-solar planets 



with regard to their tidal evolution (jRasio et al.lll996l : ISasselovll2003l : iTerquem 
2007h . 



1998 



Qgilvie fc Lin 



Finally, we note that the stellar flux incident on WASP-19b at the substellar p oint is 3.64 x 



10 9 er gs cm s placin g it in the class o f highly irradited plan ets like OGLE-TR-56b ( Konacki et al 



2003 



Torres et al. 



2004 ) and WASP-lb (jCameron et al.ll2007l ). Therefore, we expect the planet to 



have large secondary eclipse depths in the mid-IR, evidence of an at mospheric temperat ure in- 
version, molecular emission features, and a large day/night contrast (jFortnev et al.1 120081 ) . The 
existence of a hot stratosphere can be tested using secondary eclipse measurements that are cur- 
rently being obtained with Spitzer. Ec lipse measurements in the near-IR and optical z-band are also 
possible given current technology (e.g. lGillon et al.ll2009bl ). Furthermore, the density of WASP-19b 
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is half that of Jupiter's (p p = 0.51pj), so it is slightly bloated for its mass, but not extremely so. The 
high irradation, increased metallicity of the host star, and dis s ipation of tidal energy are possible 



facto rs causing the enhanced radius (IBodenheimer et al.l 120031 ; iFortnev et al.l 120071 ; iBurrows et al 
2007h . 



In summary, WASP-19b is the shortest period transiting planet yet detected. It has a mass, 
M p i=1.15 Mj and radius, -R p i=1.31 Rj. The planet orbits a main sequence G-dwarf with a slightly 
super-solar metallicity and a rotation period of P ro t = 10.5 ± 0.2 days. It is likely WASP-19b 
has been spiraling into its host star over its lifetime and has spun up the star in the process. A 
more precise age determination of WASP- 19 will allow us to confirm this and to place stronger 
constraints on the star's tidal quality factor, Q' s . 



The SuperWASP Consortium consists of astronomers primarily from the Queen's University 
Belfast, St Andrews, Keele, Leicester, The Open University, Isaac Newton Group La Palma and 
Instituto de Astrofi'sica de Canarias. The SuperWASP Cameras were constructed and operated 
with funds made available from Consortium Universities and the UK's Science and Technology 
Facilities Council. 
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Fig. 1. — WASP-South discovery photometry of WASP-19. The data are phase-folded with the 
ephemeris given in Table El 
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Fig. 2. — FTS z-band photometry of the WASP-19 transit. The data are converted to phase using 
the ephemer is given in Table [ 3l Oy erplotted is the best fitting model transit light curve us ing the 
forma lism of iMandel fc Agoll (|2002l ) applying the 4th-order limb darkening coefficients from IClaret 
( 2004h . 
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Phase 

Fig. 3. — Radial velocity curve of WASP- 19 phase-folded with the ephemeris given in Table |3l 
Overplotted is the best fitting model curve obtained from a combined analysis of the photometric 
and spectroscopic data. 
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Fig. 4. — Line bisector span versus radial velocity. The uncertainties on the line bisector values 
are double the values on the radial velocity measurement. There is no correlation between the line 
bisector span measurements and radial velocity which rules out star spot variations or a blended 
eclipsing binary as the cause for the velocity variations. 
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Fig. 5. — Modifed Hertzprung- Russell diagram c omparing the stellar density and temperature of 
WASP- 19 to theoretical stellar evolution tracks bv lGirardi et al.1 ((2000) interpolated at a metallicity 
of [M/H]=+0.1. The solid circle shows the result when the eccentricity is zero, and the asteriks 
shows the value for the stellar density if the eccentricity is a free floating parameter. The mass 
tracks are labelled and the isochrones are 0.1 (solid), 1 (dashed), 5 (dot-dashed), and 10 (dotted) 
Gyr. 
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Fig. 6. — WASP-South light curve data from 2007 phase-folded on the rotation period detected in 
the sine fitting, P rot = 10.5 days. 
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Fig. 7. — Top: Periodogram, Ax 2 /Xbest versus frequency, resulting from fitting a sine wave to the 
2007 WASP-South light curve. The peak in the periodogram has a FAP = 1 x 10 -7 and occurs at 
Prat = 10.5 days indicating the rotation period of the star. Bottom: Periodogram resulting from 
fitting a sine wave to the shuffled 2007 light curve in which the periodic signal was destroyed, but 
the red ([Pont et al.l 120061 ) and white noise were preserved. We adopted a FAP = 0.95 for the 
highest peak. 
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Fig. 8. — Orbital separation versus age of WASP-19b when evolved forward from the current state to 
the point when the planet overflows its Roche lobe. Due to exchange of angular momentum through 
tides, the planet loses orbital angular momentum and spirals into the star. The different lines which 
are labelled correspond to different values for the stellar quality factor, Q' s . The remaining lifetime 
of the planet ranges from 4 Myr to 4 Gyr depending on the value of Q' s . 
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Fig. 9. — Stellar rotation versus age of WASP-19b when evolved forward accounting for tidal 
evolution. The star gains rotational angular momentum and spins up as the planet spirals into 
the star, losing angular momentum. The different lines are labelled for different values of Q' s . For 
Q' s = 10 9 , the tidal evolution is slow enough that magnetic breaking due to stellar winds dominates 
the angular momentum of the star causing it to slow down prior to the eventual spiral in of the 
planet after ~ 4 Gyr. 



